FERMILAB-Pub-00/147-T 
OSU-HEP-00-03 
CTP-TAMU-00-20 
UMD-PP-00-084 



Enhanced Electric Dipole Moment of the Muon 
in the Presence of Large Neutrino Mixing 



K.S. Babu 1 2 , B. Dutta 3 and R.N. Mohapatra 4 

1 Theory Group, Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 
2 Department of Physics, Oklahoma State University, Stillwater, OK 74078, USA* 
3 Center for Theoretical Physics, Department of Physics, Texas A & M University, 

College Station, TX 77843, USA 
4 Department of Physics, University of Maryland, College Park, MD 20742, USA 



Abstract 

The electric dipole moment of the muon (dZ) is evaluated in super symmetric 
models with nonzero neutrino masses and large neutrino mixing arising from 
the seesaw mechanism. It is found that if the seesaw mechanism is embedded 
in the framework of a left-right symmetric gauge structure, the interactions 
responsible for the right-handed neutrino Majorana masses lead to an en- 
hancement in d e „. We find d e „ as large as 5 x 10 -23 ecm with a correlated value 
of (g — 2)^ ~ 13 x 10~ 10 , even for low values of tan/3. This should provide 
a strong motivation for improving the edm of the muon to the level of 10 - 
ecm as has recently been proposed. 
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I. INTRODUCTION 



It has long been recognized that electric dipole moments (edm) of fermions can provide 
a unique window to probe into the nature of the forces that are responsible for CP violation 
H[. Experimental limits on the edm of neutron have reached the impressive level of 6 x 10 -26 
ecm and have already helped constrain and sometimes exclude theoretical models of CP 
violation. Currently efforts are under way to improve this limit by at least two orders of 
magnitude [|3] , which will no doubt have very important implications for physics beyond the 
standard model. Electric dipole moment of the electron has severely been constrained by 
atomic measurements in Cs {d e e < 10~ 26 ) and Tl [d e e < 4.3 x 1CT 27 ecm) The limits on 
the muon edm on the other hand are much weaker, the present limit derived from the CERN 
(g — 2) experiment |5| is d e < 1.1 x 10~ 18 ecm. There has been a recent proposal to improve 
this limit on d^ to the level of 10~ 24 ecm ||. In this paper we will argue that there is a strong 
motivation for this proposed improvement, related to the observation of neutrino masses and 
oscillations. We will show that a natural understanding of small neutrino masses with large 
oscillation angles in the framework of the seesaw mechanism will lead to an enhancement 
of d% to values as large as 5 x 1CT 23 ecm, which is well within the reach of the proposed 
experiment. 

As for the theory of leptonic edm, in a large class of models a generic scaling law holds, 
given by d e /d\ ~ m^/m e . If such a relation is valid, even prior to any detailed calculation, 
one can infer that the present upper limit on electron edm will constrain the muon edm 
to be less than about 10~ 24 ecm. This scaling law arises due to the chiral structure of the 
edm operator, which is very similar to the operator corresponding to the fermion mass. 
To the lowest order in the light fermion Yukawa couplings, the edm becomes proportional 
linearly to the fermion mass. In specific models, it may so happen that other constraints 
put the electron edm itself at a much lower value; e.g., the standard model prediction for 
the electron edm is ~ 10~ 41 ecm 0. The scaling law then suggests that the corresponding 
value for the muon edm would be at the level of 1CT 39 ecm, which is beyond the reach of 
any conceivable experiment. In multi-Higgs doublet extensions of the standard model, the 
dominant contribution to the leptonic edm arises from a two-loop diagram involving 7 — V- 
Higgs vertex, where V = Z, W || . Since such a vertex is flavor universal, when converted 
to the fermion edm, the above-mentioned scaling law will hold. Recently an extended Higgs 
model H has been analyzed, where it has been shown that for large values of the parameter 
tan/3 (ratio of the two Higgs vacuum expectation values), the one-loop diagram that scales 
as m^A 2 , where A^ is the muon Yukawa coupling, can compete with the two-loop diagram 
||, leading to order one violation of the scaling law. 

In the supersymmetric extension of the standard model (MSSM), under the usual as- 
sumptions about supersymmetry breaking terms, i.e., universality of scalar mass terms and 
proportionality of the trilinear A terms with the corresponding Yukawa couplings, a similar 
scaling law would hold. A leading contribution to leptonic edm in such models is the one- 
loop diagram involving the bino virtual state and a complex An term. The assumption of 
proportionality of A terms then implies that the above mentioned scaling relation remains. 
A similar remark holds when the chargino diagram is considered, with a complex fi term, 
again due to the universality of the CP violating parameter. (For a discussion of edm of 
electrons in MSSM and SUGRA models, see ref. [p!0|-[l3[|.) Evaluation of these bino and 
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chargino diagrams leads to a value for the muon edm of about 8 x 10 -25 ecm, once the 
upper limit on electron and neutron edm are satisfied. The expected reach of a proposed 
BNL experiment for the muon edm is 10~ 24 ecm, which is somewhat above the largest value 
allowed within the MSSM. 

Recent experimental evidence for neutrino masses, especially from the SuperKamiokande 



atmospheric neutrino data ||14|| , suggests that the MSSM must be extended to account for 
it. A natural place for small neutrino masses is the left-right symmetric extension of the 
standard model jnj. We have recently advocated a simple supersymmetric realization of 



left-right symmetry (SUSYLR) which accommodates the neutrino masses via the seesaw 



mechanism | IB| . Our proposal is simply to embed the MSSM into a left-right symmetric 
gauge structure at a high scale v r ~ 10 11 — 10 15 GeV. The effective MSSM that emerges from 
this model at scales below the left-right symmetry breaking scale, v r, is a constrained MSSM 
with far fewer number of phases. In particular, it has a built-in solution to the SUSY CP 
problem piTyTBf. Owing to the constraints of parity symmetry, the Yukawa coupling matrices 



and the trilinear A matrices become hermitian in this model. Similarly, the /z term, the soft 
B[i parameter, and the gluino mass parameters all become real, eliminating potentially 
excessive CP violation from the MSSM. Furthermore, i?-Parity arises automatically in this 
model as part of the gauge symmetry, since the gauge structure involves B — L symmetry. 

In this paper we wish to investigate the CP violating muon edm and (g — 2) M in this 
class of models. We will show that the interactions responsible for the Majorana masses of 
the right-handed neutrinos will lead to an enhancement of d^. We find d e ^ as large as 5 x 10~ 23 
ecm and (g — 2) M as much as 13 x 1CT 10 . These values arise even for small tan/3 ~ 3. Our 
main effect arises through the renormalization group extrapolation from the Planck scale 



to the left-right scale vr In this interval the Yukawa couplings of the vr fields which 
induce their Majorana masses, as well as the associated trilinear A terms, will affect the soft 
supersymmetry breaking parameters of the effective MSSM, leading to the enhancement of 
d^. Sine the Majorana Yukawa couplings do not obey e — fi universality, the scaling law 
^il/^e = m n/ m e is not obeyed by these new diagrams. 

For concreteness, we will work within the framework of a minimal version of the high 
scale SUSYLR (or 5*0(10)) model. It is minimal in the sense that we have only one multiplet 
of Higgs field that gives rise to the usual Dirac fermion masses, i.e., one left-right bidoublet 
$ (10 in the case of SO(10)). With one such multiplet, only one Yukawa coupling matrix is 
allowed in the quark sector, leading to the proportionality of the up and the down Yukawa 
coupling matrices 0,^,^T|] . We call this up-down unification. It has the consequence that 
all the flavor mixings vanish at the tree level. We have shown that acceptable values of the 
mixing angles can arise from the one-loop diagrams involving the gluino (and the chargino), 
proportional to the flavor structure of the trilinear A terms. This considerably restricts the 
flavor and CP violating interactions in the model and makes it very predictive. The model 
has been shown to lead to a consistent picture of Kaon CP violation including e and e' and 
it predicts neutron edm at the level of 10~ 27 ecm. The leptonic sector of the model was 
investigated in Ref. f2"Uf , we shall work within that framework to calculate the edm and 
(g — 2) of the muon. We have verified that going to non-minimal models, e.g., by employing 
more than one bidoublet Higgs field, does not affect our results by much in the leptonic 
sector. 
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II. BRIEF OVERVIEW OF THE MODEL 



Let us briefly review the salient features of the minimal SUSYLR model. The electroweak 
gauge group of the model is SU{2) L x SU{2) R x U{1)b~l with the standard assignment of 
quarks and leptons - left-handed quarks and leptons (Q, L) transform as doublets of SU(2) L , 
while the right-handed ones (Q c , L c ) are doublets of SU(2)r. The Dirac masses of fermions 
arise through their Yukawa couplings to a Higgs bidoublet $(2, 2, 0). The SU (2) R xU(l) r-l 
symmetry is broken to U(l)y by B — L = 2 triplet scalar fields, the left triplet A and right 
triplet A c (accompanied by A and A c fields, their conjugates to cancel anomalies). These 
fields also couple to the leptons and are responsible for inducing large Majorana masses for 
the vr. An alternative is to use B — L = 1 doublets x (left) and \ c (right) along with \ 
and x c instead of the A fields. Here we shall adopt the B — L = 2 triplet option, which 
allows direct couplings to the leptons and which conserve i?-Parity automatically. Let us 
write down the gauge invariant matter part of the superpotential involving these fields: 

W = Y q Q T T 2 $T 2 Q c + Y^ r r 2 $r 2 L c 

+ (fL T zr 2 AL + i c L cT ir 2 A c L c ) . (1) 

Under left-right parity, Q <-> Q C *,L <-> L c *,$ <->• $+, A <-> A c *, along with W S u{2) L 
Wsu(2) R i Wb-l *-* W B _ L and 6^6. Here the transformations apply to the respective 
superfileds. As a consequence, Y q = Yj, Y/ = Yj, and f = f * in Eq. (1). Furthermore, the 
trilinear A q and A\ terms will be hermitian, gluino mass term will be real, and the supersym- 
metric mass term for $ (the /i-term) as well as the supersymmetry breaking Bfi term will be 
real. Departures from these boundary conditions below vr due to the renormalization group 
extrapolation is small. The model thus provides a natural resolution to the supersymmetric 
CP problem. 

Below Vr, the effective theory is the MSSM with the H u and H d Higgs multiplets. These 
are contained in the bidoublet $ of the SUSYLR model, but in general they can also reside 
partially in other multiplets having identical quantum numbers under the MSSM symmetry. 
Allowing for such a possibility, the single coupling matrix Y g of Eq. (1) describes the flavor 
mixing in the MSSM in both the up and the down sectors leading to the relations 

Y u = 7 Y d , Y, = 7 Y v d, (2) 

which we call up-down unification. Here 7 is a parameter characterizing how much of H u 
and Hd of MSSM are in the bidoublet $. The case of H U ( i entirely in $ will correspond to 
7 = 1 and tan/3 = m t /mb. At first sight the first of the relations in Eq.(2) might appear 
phenomenologically disastrous since it leads to vanishing quark mixings and unacceptable 



quark mass ratios. We showed in Ref. |L8| that including the one-loop diagrams involving 
the gluino and the chargino and allowing for a flavor structure for the A terms, there exists 
a large range of parameters (though not the entire range possible in the usual MSSM) where 
correct quark mixings as well as masses can be obtained consistent with flavor changing 



constraints. In Ref. ||18|| , we explored the parameter space that allowed for arbitrary squark 
masses and mixings as well as arbitrary form for the supersymmetry breaking A matrix. We 
found a class of solutions for large tan j3 ~ 35 — 40 (7 =1), and for small tan j3 ~ 4 where all 
quark masses mixings and CP violating phenomena could be explained. The smaller value 
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tan/3 requires larger values of 7, since 7 tan/? = m t /mb is fixed. In this paper, we use small 
tan (3 scenarios which is less constrained. 

Since the parameter 7 plays a crucial role in determining the value of tan f3, let us explain 
its origin in an explicit high scale model. We will also show how the solution to the SUSY 
CP problem can be maintained even for the case of small tan f3. 7 arises from the mixing 
of the bidoublet $ with other weak doublets in the high scale theory. We assume that only 
one pair of doublets, H u and Hd of MSSM, remain light below vr. A concrete example 
which also maintains automatic i?-Parity of the left-right model involves the addition of 
the following new fields: p(2, 2, 2) + p(2, 2, -2) and Q L (3, 1, 0) + Q R (1, 3, 0). They lead to 
the following new terms in the superpotential: 

Wnew = m(aa + a c a c ) + /^$ 2 + fi p pp + im(n 2 L + n%) 

+ Ai [Tr(pA c $) + Tr(pA$)] + A 2 [Tr(pA c $) + Tr(pA$)" 

+ A 3 Tr(A c fi fl A c + AQ L A) + X 4 Tr(pft R p + pQ L p) . (3) 

The coupling and the mass parameters in Eq. (3) are guaranteed to be real by parity 
symmetry, P, defined earlier in combination with the charge conjugation symmetry C under 
which all superfields (except p and $) transform as ^ — > \l/ c , where \I/ stands for a relevant 
superfield in the theory; the Wl Wr and B — > —B. The fields p and $ transform as 
follows: p — > r 2 p T T 2 and $ — > r 2 $ T r 2 . We will assume that the supersymmetry breaking 
terms respect only P and not C . 

It can be shown (see e.g., Ref. ||26|| ) that this model has a ground state where {Qr) ~ 
(A c ) = (A c ^> ~ vr and (Ql) = (A) = (Ay = 0. The p superfield contains an H u -\ike 
MSSM doublet and p contains an if^-like one. Once the right handed gauge symmetry is 
broken by A c vev, the doublets in $ and those in p and p mix via a matrix, which is given 

by W mass = ( p u $ u ) M doub i e t ( ^ ) , where 



Mdoubiet being an asymmetric matrix leads to light eigenstates 
given by H u = cos^i$ u + sin^ip^ and Hd = cos 9 2 &d + sm9 2 pd- Here 6\ is the 
Pu — mixing angle, which is unrelated (due to the asymmetry of the matrix) to 9 2 , the 
Pd — &d mixing angle. This gives 7 = ^ * , which can take any arbitrary value. 

We note that due to the combination of P and softly broken C symmetry, all dimension 
four couplings are real, leading to a solution to the SUSY CP problem. To see this, note that 
due to these symmetries, all entries in the mass matrix of Eq. (4) are real, so that the effective 
p term of MSSM stays real. (With parity symmetry alone, the Ai j2 couplings in Eq. (4) could 
be complex, which would make the effective p term of the MSSM complex.) Furthermore, 
since only the dimension 3 and 2 terms of the SUSY breaking Lagrangian are assumed to 
respect P, but not C, such a scenario is completely stable under renormalization. (This 
scheme is distinct from scenarios where CP symmetry is imposed on the MSSM Lagrangian 
at a high scale to solve the SUSY CP problem El. Since the gauge structure of MSSM 



does not have parity symmetry, the phases of the soft SUSY breaking terms will have to be 
small in that case.) 
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Unlike the large tan/? case (corresponding to 7 = 1), we are finding that CP violation 
in the quark sector has to arise from soft terms. We have analyzed this possibility in Ref. 
I8| and shown its consistency. We are pursuing this possibility further |23j. An immediate 



outcome of this scenario for hadronic CP violation is that although there is KM type CP 
violation, generically it tends to be sub-leading to SUSY CP violation. 

In the absence of the Qr field in Eq. (3), the doubly charged field A c++ in A c (as well 
as A c in A c ) will remain massless - it will pick up mass only of order the weak scale, or 
of order v R / M string , if non-renormalizable operators are included. Inclusion of Qr (and its 
left-handed partner fl^) lifts the mass of A c++ to the scale vr [J2b|. We will analyze two 



cases, one with the inclusion of Vt L ^ R fields, and one without. In the latter case, we will take 
the mass of A c++ to be ~ v\j 'M string . 



III. LEPTONIC CP VIOLATION AND MUON EDM 

To discuss CP violation in the lepton sector, we need to specify the leptonic superpoten- 
tial Wt and the most general soft breaking Lagrangian, C e soft , in the lepton superpartners. 
The leptonic Wt is given in Eq. (1), C e soft is given by: 

- £f oft = m 2 LL DL + m 2 RR L ct L c + [A^L? + A f (LLA + L C L C A C ) + H.c] (5) 

To generate a nonvanishing muon edm, one needs a complex valued (Ai) 2 2 and/or complex 
soft mass-squared terms. But above the scale where the parity symmetry is valid, Ai is 
hermitian and therefore its diagonal elements are all real. This element can however be 
complex due to radiative corrections below the parity breaking scale. There are two ways 
this can happen: (i) if only parity symmetry is broken but gauge symmetry SU{2)i J x 
SU{2)r x U{1)b-l is unbroken at the string scale by introduction of parity odd singlets 
||24|| ; (ii) if both parity and the left-right gauge symmetry are broken, but some remnant 
of the f and Af couplings remain below the vr scale. This has been shown to happen in 



supersymmetric left-right models with minimal field content [p5 |. In the explicit version 
described in Sec. II, if the flL,R fields are absent, the A c++ field from A c will have a mass 
of order v R /M stTing ~ 10 12 GeV. So between M string and M&c++, the effects of f and Af 
couplings will be felt, and (Ai) 22 can become complex. This will also induce flavor violating 
complex soft mass-squared terms proportional to AfAL, even if we start with diagonal soft 
masses at M string . 

In case (i), the way (Ai)ij become complex is as follows. Below the D-parity (discrete 
parity) breaking scale, M string , only A c, s (and not A's) contribute to renormalization group 
equations (RGE) describing the evolution of (Ai)^ since the A's acquire masses of order 
String- The RGE are given in the Appendix for this case. We have, from Eq. (25) of 
Appendix, 

^I^f'Y.^ffH). (6) 

The first term on the RHS of Eq. (6) will introduce phase in A\. Note that Af is not 
constrained to be hermitian at the string scale by parity symmetry (unlike Ai, which must 
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be hermitian at M str i n g). We will allow for complex entries in the 2 — 3 block of Af in our 
analysis. 

Below the D-parity breaking scale, the soft mass parameters m| L and m 2 RR will evolve 
differently. In particular, m RR will feel the effects of f and Af couplings. In order to explain 
the large oscillation angle needed for the atmospheric neutrino data, we will find that /2 3 
is not much smaller than / 33 . Thus (Af) 2 3 is not much smaller than (A/) 33 . Consequently, 
(m RR )23 will become large and complex. This is the main source of the enhanced edm of 
the muon in the model. This qualitative feature becomes more transparent if we examine 
the RGE for m RR (see Eq. (27) of Appendix). It has the form: 



dm 2 



RR oc 



dt 167T 

It is clear from Eq.(7) how (111^)23 becomes large and complex. 

The dominant contribution to the edm of muon arises from a diagram which has right 
and left-handed muon in the external legs and a lighter stau inside the loop. It utilizes the 
above-mentioned 2 — 3 mixing which is large and complex. For example, the diagram can 
have hl — tr and f L — fi R vertices along with the stau mass flip inside the loop or it can 
involve just the — f R and r R — ji R vertices. It might be suspected that similar diagrams 
will also induce large edm for the electron. However, in this model, since /i 3 and /12 are 
much smaller, such contributions are negligible. Essentially, we have a scenario where e — /i 
flavor symmetry is broken by a large amount by the f and Af terms. As a result the scaling 
law alluded to in the introduction does not hold. If we assume, as we do in our analysis, 
the existence of phases only in the 2 — 3 block of the Af matrix, or if Af has negligible 
entries in its first row and column, no appreciable edm for the electron gets induced due to 
mixing effect. Below v R , we have only the MSSM field content. Due to the new f couplings 
above the v R scale, the fi mass is lower than usual SUGRA model for the same values of the 
parameter space (i.e., m , m^, A Q , tan/?). This is why the diagram involving the f tends 
to dominate in d e ,,. 

In case (ii), we use the fact that in the minimal SUSYLR model (without Ql,r), A c++ 
and A c remain below the v R scale; therefore their couplings to the charged fermions via 

RGE's lead to imaginary parts in (A/)22 by an amount '—^r—- Again the soft masses 
become complex in the same fashion as in case (i). These fields get decoupled at somewhat 
lower scale ~ 10 11 GeV, below which the spectrum is that of MSSM. 

IV. RESULTS 

Let us first discuss the neutrino mass fits in this model. We start with a basis where the 
charged leptons masses are diagonal and Dirac neutrino masses are given by 

M u d = 7 tan (3 Mi, (8) 

where Mi = Diag(m e ,m[t,m T ). The light Majorana neutrino mass matrix is then given by: 

M u = 1 ^l Ml r 1 M h (9) 
vr 
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where f is the right-handed Majorana Yukawa coupling matrix. 

In our fit, we first use the small angle MSW oscillations solution for the solar neutrino 
deficit with AmJ, ~ (0.3 - 1) x 10~ 5 eV 2 and 2 x 10" 3 < sin 2 2# eM < 2 x 10" 3 . We also use 
the Up — > v T oscillation scenario to explain the observed deficit in the flux of muon neutrinos 



from the atmosphere |HJ]. The mass splitting is taken to be Am 2 ~ (0.1 — 1) x 10 eV 



and the oscillation angle to be sin 26 ^ T ~ 0.8 — 1. 

For tan/3 = 3, we find a good fit to the solar and atmospheric neutrino data by choosing 
f at M string to be 

/-l.OOxlO" 4 8.8 x 10~ 4 -2.2xl0~ 5 \ 
f= 8.8 x 10~ 4 -1.3 x 10~ 2 1.03 x HT 1 . (10) 
\-2.4xl0~ 5 1.03 x lO" 1 -1.59 / 

The resulting neutrino masses at vr = 10 15 3 GeV are: (6.27 x 10~ 6 , 2.5 x 10 -3 , 5.2 x 10 -2 ) 
eV. The leptonic mixing matrix is given by: 



U 



I -0.99 

3.1 x 10~ 2 0.74 -0.67 | . (11) 

V -2.9 x 10 



-2 




U21 is the mixing angle relevant for solar neutrino oscillations. (Our notation is such that 
UM V U T = j\/f diagonal ^ This choice leads to a simultaneous explanation of the solar and 
atmospheric neutrino anomalies. Note that we have taken all Yukawa couplings to be real, 
consistent with our assumption that C and P symmetry are respected by d = 4 terms. 

It is possible to fit the large angle oscillations solution to satisfy the solar neutrino deficit. 
In that case we take f matrix is at vr ~ 10 15,6 GeV to be 



/-1.77xl0~ 7 -1.42 x 10" 6 
f = I -1.42 x 10~ 6 -3.9 x 10~ 3 

-6.4 x 10~ 2 -1.28 





-6.4 x 10~ 2 I . (12) 



With these values the neutrino masses are (1.7 x 10 3 , 2.0 x 10 3 , 3.4 x 10 2 ) eV and the 
corresponding leptonic mixing matrix is: 

/ -0.89 -3.3 x 10- 1 -4.0x10-^ 
U= -4.5 xlO- 1 0.63 0.63 . (13) 

\-2.8xl0- 2 0.72 -0.69 / 

We use the one-loop Yukawa and two-loop gauge RGE to extrapolate all parameters 
between the string scale and the vr scale. Since the new couplings f affect the RGE for the 
leptonic Yukawas, one needs to make sure that the charged lepton masses come out to be 
correct at the weak scale. For simplicity we choose a universal scenario, i.e., all the scalar 
masses are given by a common mass parameter m at the string scale. We also assume a 
common trilinear mass Aq(xYi) for all generations. For Af we use a structure similar to f. 
But we do not impose Af oc f. We demand electroweak symmetry to be broken radiatively. 
In case (i), where parity is broken at M string , A fields get decoupled and only the A c fields 
contribute to the RGE for soft masses. Consequently the renormalized right handed slepton 
masses get lowered due to the presence of the new couplings f . Furthermore, Ai will pick 
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up off-diagonal elements and will lose its hermitian structure through renormalization. The 
A c fields get decoupled at the left-right breaking scale Vr, below which we use the RGE 
corresponding to the MSSM degrees. 

The EDM for a spin 1/2 fermion is given by the following effective Lagrangian: 

L f = —d^a^F^ . (14) 

In this model, we have only the neutralino-slepton loop contribution to the edm of muon. 
This contribution is given as |12| |: 

M Airsin z 6 w ~{~{ m k m k 

where 77^ = [-v^(tan6' H /(Q M -r3 M )X H + r3 M X2i)r^ 2fc +x^X 3i r|j 2fc ](v^tan6' w <5 At Xi i r i? 2fc- 
x M X 3 jr^ 2 fc) and x M = ^^p*^ ^ • -X" diagonalizes the neutralino mass matrix, X T M x oX = 
diag(m x o,m x o,m x o,m x o 4 ). Here Tl,r are 6x3 matrices given by qi = ^l,rQ and B(r) = 
w ± w (l + r + 2r\nr/(l + r)). 

We first analyzed the case where Af and f are proportional. It still allows for an overall 
phase in Af, consistent with P invariance. In this case d^ is highly suppressed, d^ < 1CT 26 
ecm. The reason is that with only one matrix structure f, when the effective (Ai) 2 2 is 
computed in the original gauge bases, it will remain real. Small contribution will arise in 
the mixed \x — r EDM operator, which can lead to a small value of d^ since the physical 
fj, is a linear combination of the two states. However, this fi — r mixing turns out to be 
small. As soon as the proportionality Af oc f is relaxed, d^ becomes much larger. We 
have analyzed the case where Af and f are non-proportional, but the magnitudes \Af\ij are 
proportional to |f|^. We allow phases of order 1 in the (23), (32) and (33) elements of Af 
matrix, while keeping real. In this case we find the maximum muon edm to be 7 x 1CT 25 
ecm. When this assumption of proportionality of the magnitudes is relaxed, even larger 
value of d^ results. We give an explicit example for this case below. It should be mentioned 
that large values of Af reduces stau mass while it increases d^. So in exploring regions of 
large d e , we need to consider the experimental limits on stau. In our calculation we take 
the lightest stau mass (fi) to be > 80 GeV (which is above the current experimental limit 
of 70 GeV []27| at y/s = 202 GeV). In Fig. 1, we exhibit the case which has small angle 



oscillation solution. The large angle solution, however, does not show any difference. In Fig. 
1 we plot the muon edm parameter = £ogio[ lxlo -23 ecm ] for case (i) for tan/? = 3. This 
corresponds to D-parity broken at the string scale, but left-right gauge symmetry broken 
at vr ~ 10 15 ' 3 GeV. At the string scale (taken to be 10 17 GeV), we have assumed (in GeV 
units throughout) 

/ -2xl0- 3 1 x 10~ 2 \ 

A f = 1.0 x 10~ 2 -1 x 10 2 e^ /2 4.7 x 10 2 e i7r/2 . (16) 
V 4.7 x 10 2 e i7r / 2 3.3 x 10 2 e-^/ 2 ; 

We put Aq = —120 GeV (with A\ = A Yi). The solid line in Fig. 1 is drawn for mo=160 
GeV. The extreme left corner of the curve corresponds to lighter stau mass (fi)=82 GeV. At 
the same spot in the parameter space, the lightest chargino (xt) and the lightest neutralino 
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masses (xi) are 106 GeV and 52 GeV respectively. We can see that the muon edm can be 
as large ~ 3 x 1CT 23 ecm in this case. The dotted line is drawn for m =170 GeV for the 
same set of input values. 

In Fig. 2 we plot the muon edm parameter for case (ii) with tan/? = 3 and mo = 160 
GeV. This case corresponds to A c++ surviving below vr. We assume the scale at which it 
decouples to be 10 12 GeV. We have used the universal scenario for the slepton masses and 
have used the same f matrix as before. At the string scale, we take (in GeV) 

/-2xl0~ 3 1 x 10~ 2 \ 

Af = 1.0 x 10" 2 -1 x lOV*/ 2 3.0 x 10 2 e i7r / 2 . (17) 
V 3.0 x 10 2 e i7r / 2 1.1 x 10 2 e~ i7r / 2 y 

We take Aq = GeV. The extreme left corner of the curve in Fig. 2 corresponds to lighter 
stau mass (fx) mass of 80 GeV. At the same spot, as before, the Xi an d the Xi masses (x?) 
are 106 GeV and 52 GeV respectively. As can be seen from the figure, large values of d e are 
possible, as large as 5 x 10 -23 ecm. 

We have assumed non-proportionality of Af and f in the preceding two examples. We 
will argue that this is not unnatural. First of all, there are no strong experimental hints that 
suggest proportionality of the two (unlike the case of A[ and Yj). Second, we have proposed 
recently a model based on horizontal gauge symmetry which allows for all parameters of 
the soft breaking sector to be arbitrary, subject only to the constraints of the horizontal 



symmetry H [|28]. The symmetry H was taken to be SU(2)h x U(1)h, with the first two 
generations of fermions falling into SU(2)h doublets and the thrid generation into singlets. 
The first two generations have U(1)h charges of —1, while the third generation is neutral. 
H is spontaneously broken by a pair of doublet [0(+l), 0(— 1)] and singlet [x(+l), x(— 1)] 
scalars fields whose vev's are below the string scale. We denote = (<p) /M string , e x = 
(x) /M s t r ing with ~ 1/7, e x ~ 1/25. The effective Yukawa couplings involving the light 
two generations will be proportional to powers of and e x . The U{1)h also alleviates 
potential problems with D-terms associated with horizontal symmetries. 

Within the SU(2)h x (U1)h model, it is not necessary to assume universality of scalar 
masses or proportionality of A terms and the Yukawa couplings. For the first two generations, 
the scalar masses will be approximately equal, owing to the non-Abelian sector of the 
horizontal symmetry. With the horizontal charge assignment given above, we can write 
down the most general if-symmetric Yukawa couplings, soft mass terms and A terms. 
Since the A terms become hierarchical, all FCNC constraints can be satisfied, even without 
proportionality assumption |28| . 



We will now present an example for the muon edm within this horizontal symmetric 



framework. We will embed the model of Ref. [28] into left-right symmetry at a high scale. 



Unlike in Ref. ||28|| , all the CKM mixing will vanish at tree-level now. In a basis where the 
Yukawa couplings are diagonalized, the Majorana neutrino coupling can be written in the 
following hierarchical form: 

/n^/e^ fi2 e x fiz e x / e <l>\ 

fi2t x / 22 4 fwH • ( 18 ) 
,/l3 e x / e ^ /23 e ^33 / 

The bilinear soft mass matrix and the A matrix are given as: 
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ni 



RR 



rrir 



rrin 



m 2 {x l2 )e x 
m§(xi 3 )4/e,£ mg(x 23 )e0 



m 2 (x 12 )e x m 2 (x 13 )e 2 x / e/ 



rn 



:« 



A = An 



(2/21) e x 
V(2/3i)4/^ 



(2/22)4 
(?/32)e<^ 



(yis)4/ e 
(2/23)^ 
2/33 



(19) 



We also have m 2 LL = m 2 RR . This structure for A hold for both A\ and A/ (as well as for A q ). 
At Muring we will take A\ to be hermitian. In order to fit the experimental values of quark 
and lepton masses we choose = 1/7 and e x = 1/25. In this new scenario, the muon edm 
can be enhanced to 5 x 10 -23 . We have taken soft masses for all the Higgs fields to be 85 
GeV. In Fig. 3, we exhibit the results for d e ^ for one such example. To generate this plot, 
the input values we have used at the string scale are as follows: 



/ 1 



rn 



RR 



85 2 



2 fc X 



2 X 
1 



1.8 



; A, = 30 



<l*4 

2 t X 



±e 2 

3 x 



2 fc X 
if 2 

—in/3. 



ie 2 

3 C X 

-4 



(20) 



f = 



/ -1.01 x HT 4 
9.0 x 10~ 4 
-1.4 x 10~ 3 



V 



A f = 500 



3 e x 



v 7 X 



9.0 x 10~ 4 -1.4xl0-' 3 ' 
-1.2 x 10~ 2 1.04 xlO" 1 
1.04 x 10" 1 -1.59 

3 fc X 

-3e^/ 2 e 2 

4e i7r /% 0.6e~ i7r / 2 




(21) 



Note that we have allowed for all coefficients to be order one, consistent with the horizontal 
symmetry. (This is also true for the f matrix elements.) The (13), (31) elements of the f 
are no longer very small like our previous example because of the symmetry requirement. 

The choice of f matrix in this case corresponds to the following light neutrino masses: 
(6.27 x 10~ 6 , 2.9 x 10 _ 3 5 4.4 x 10~ 2 ) eV. The corresponding leptonic mixing matrix is: 



U 



( -0.99 4.2 x 10~ 2 

3.2 x 10~ 2 0.70 
\-2.9xl0~ 2 -0.74 




(22) 



Note that in this example v e — oscillation explains the solar neutrino data via small 
angle MSW oscillation, — v T oscillation explains the atmospheric neutrino data. We have 
found that by varying the order one couplings slightly, it is also possible to obtain a different 



oscillation is relevant for solar neutrinos, while v.. 



scenario whre v e — v 7 
with m v > m Ur explains the atmospheric neutrino data [28 
much altetered in such a scenario. 

Now we turn to the evaluation of (g — 2) of the muon. 
contribution from the chargino and neutralino diagrams 
can be found in Ref, 



oscillations 



The predictions for d e is not 



In MSSM, the (g - 2)„ gets 



2J. The relevant expressions 
In this model we have contributions from both these loops. The 
chargino contribution is somewhat bigger than the neutralino loop. We find the magnitude 



31 
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of (g - 2)„ to be (6 - 10) x 1(T 10 for the curves in Figs. 1 and 2 and (8 - 13) x 1(T 10 for 
the model with horizontal symmetry given in Fig. 3. 

As for other rare processes, the branching ratio of r — > /ry is one to two orders of 
magnitude below the present experimental limit. Since this process cnnot be made much 
smaller, it will be of great interest to improve the present limit by two orders of magnitude, 
which does not appear to be out of question. In all cases that we studied, the edm for 
electron is of order 10~ 28 ecm. As for fi — > cy, it is three to four orders of magnitude smaller 
than current limits for cases (i) and (ii), and one order of magnitude smaller than current 
limits in the case of horizontal symmetry. 

In conclusion, we have shown that in supersymmetric extensions of the stanadard model 
that accommodates small neutrino masses via the seesaw mechanism, there is an enhance- 
ment of the muon electric dipole moment. Interactions responsible for the generation of 
Majoran masses for the right-handed neutrinos are responsible for this enhancemnt through 
renormalization group effects. We have found values of d e ^ as large as 5 x 10~ 23 ecm. Our 
finding should provide a strong motivation to improve the limit of d e to the level of 1CT 24 
ecm, as has recently been proposed. Probing d e ^ at this level could reveal the underlying 
structure responsible for CP violation as well as for the generation of neutrino masses. 
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VI. APPENDIX 



In this Appendix we give the renormalization group equations appropriate for the mo- 
mentum range between M string and Vr for the case where parity is broken at M string . 

ft = vh [ ~ A7T{7aR + l aB ~ L)1 + 2f ft + 4Y ' Y ? + Tr ( f ft )] f ' ( 23 ) 
IF = jt^ [ ^ 3aR + 3aL + \ aB - L)1 + 2f ft + 4YY ? + Tr (3Y 9 Y ? t + Y,Y,t)]Y,, (24) 

ft = TfL* [~ 47r ( 3ai? + 3ai + A (25) 

I 87r(3aflM fl + 3a L M L + -a B - L M B ^ L )Yi 

+ ^YfY, + 8Y,Y,U, + 2f ft A + 4^ftY, 
+ 2Tr(AYt)Y + Tr(YYt)^ + 6Tr(A 9 Yj)Y + 3Tr(Y,Yt)Y i ], 
dAf 1 9 9 

"IT = T7T^[- 47r ( 7a « + o®B- L )A f + 87r(7a fl M fl + -a B _ L M B _ L )f (26) 

+ 8AiYji + 4YjY,U/ + 2A f tff + 4f &A f 
+ 2Tr(A / f t )f + Tr(ff t )A / ], 
dm 2 2 

= ^[-MW<*b-lM*_ l + 3a L M 2 L ) (27) 

+ i((YY/ + f tym^ + m| R (Y J Yj + f ft) + 2(Ym^Yt) 
+ 2(fm^ft + m|Y,Yj + m^f ft + A,Af + A^))]. 
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FIGURES 
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FIG. 1. k/j,(= Logxo [ 1 w 1^5 ecm ] ) is plotted against TOi/ 2 for tan/3=3 for case (i). The solid line 
is for mo = 160 GeV and the dotted line is for too = 170. The other inputs are described in the 
text. 
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FIG. 2. kfj,{= LogiQ [ 1 - 10 -^3 ecm ] ) is plotted against m-yji for tan/3=3 for case (ii). too is 160 
GeV. The other inputs are described in the text. 
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FIG. 3. k/j,(= Logio [ i x io-^3 ccm ] ) I s plotted against mi/2 for tan/?=3 for the model with hori- 
zontal symmetry, mo is 85 GeV. The other inputs are described in the text. 
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